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ABSTRACT 

Aims. After considering the effects of negative feedback on the process of star formation, we explore the relationship 
between star formation process and the associated feedback, by investigating how the mechanical feedback from su- 
pernovae(SNe) and radiative feedback from luminous objects regulate the star formation rate and therefore affect the 
cosmic reionization. 

Methods. Based on our present knowledge of the negative feedback theory and some numerical simulations, we construct 
an analytic model in the framework of the Lambda cold dark matter model. In certain parameter regions, our model 
can explain some observational results properly. 

Results. In large halos(rvir > W^K), both mechanical and radiative feedback have a similar behavior: the relative 
strength of negative feedback reduces as the redshift decreases. In contrast, in small halos (Tvir < 10'* A") that are 
thought to breed the first stars at early time, the radiative feedback gets stronger when the redshift decreases. And the 
star formation rate in these small halos depends very weakly on the star-formation efficiency. 

Conclusions. Our results show that the radiative feedback is important for the early generation stars. It can suppress the 
star formation rate considerably. But the mechanical feedback from the SNe explosions is not able to affect the early 
star formation significantly. The early star formation in small-halo objects is likely to be self-regulated. The radiative 
and mechanical feedback dominates the star formation rate of the PopII/I stars all along. The feedback from first 
generation stars is very strong and should not be neglected. However, their effects on the cosmic reionization are not 
significant, which results in a small contribution to the optical depth of Thomson scattering. 
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1. Introduction 

' In the past decades, the so-called "bottom-up" hierarchi- 
[ cal scenario for the large-scale structure formation in the 
Lambda cold dark matter (ACDM) cosmogonies has been 
getting decisive support from more and more observations 
; (e.g. from HST, WMAP, SDSS) and from high-resolution 
N-body/hydrodynamics numerical simulations. Owing to 
the complexity of the baryonic evolution in the radia- 
tive background and gravitational field, the galaxy forma- 
tion and evolution in dark matter halos and the reion- 
ization history of IGM have not been fully understood 
yet. The difficulties stem from our lack of knowledge of 
those early formed objects that are unobservable at present. 
The formation of the first objects is relevant to the cos- 
mic reionization history. Recent detecti ons of the Gunn- 
Peterson trough (|Gunn fc Peterson 19651 ) in the spectra of 
QSOs wit h z c:^ 6 indica t e less 50% neutral hydrogen at 
z ~ 6.5 ()Fan et al. 20061 : IWvithe et al. 2005D . The ongo- 
ing observations by the WMAP satellite of cosmic mi- 
crowave background (CMB) and the highest redshift QSOs 
put very tight constraints on the reionization history of 
the Universe. The WMAP five-year observation manifests 
the Thomson scattering optical depth, Te = 0.084^g Q|^g 
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(jKomatsu et al. 2008f ). which suggests that our Universe 
might be reionized during the period of redshift 9.4 < z^c < 
12.2. 

It is well-known that chemical elements heavier than 
lithium are produced exclusively through stellar nucleosyn- 
thesis. Some of the first generation stars (the so-called 
population III stars, hereafter, PopIII) die as SNe explo- 
sions, which can expel heavier elements into the intergalac- 
tic medium (IGM). When the metal elements in IGM are 
enriched to a certain threshold Z„it, the population II/I 
stars (hereafter PopII/I) will form and take the place of the 
first stars to light the universe. Thence, the SNe from PopIII 
stars determine the transition from PopIII to PopII/I. The 
existence of PopIII stars can help in explaining the metal 
enrichment from Z ~ 10"^^ — 10"^" to the lowest metallic- 
ity of PopII stars Z ~ lO"**— 10"'^, the formation of massive 
black holes, the reionization of the universe, the starting en- 
gine for t he formation of the first galaxies and the G-dwarf, 
and so on (|Ciardi fc Ferrara 20051 ). But there is still no hope 
in observing the first generation stars until the launch of the 
James Webb Space Telescope (JW ST), the s uccessor of the 
Hubble Space Telescope (HST) (|Barkana fc Locb 2001i). 
JWST is expected to find t he pair instability S Ne (PISNe) 
from massive PopIII stars (I Wise fc Abel 20051) . 
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Due to the absence of the observational data of 
very high-redshift (z > 10) objects, theoretical in- 
vestigations are mainly based on numerical simula- 
tions. Some works concentrated on the effects of 
the first-generation SNe explosions (I Yoshida et al. 2003t 
iKitavama k Yoshida 20051 : iGreif et al. 2007D . while the 
others focused o n the strong ste llar and galactic winds 
from PopIII stars (|Mevnet et al. 200 6: Ricotti et al. 2008 1. 
Based on these works, we put forward a model to describe 
the global effects from the PopIII stars at an early time, 
such as the SFR density, the IGM reionization, and so on. 

In this paper, we study the effects of the mechanical 
feedback from SNe and the radiative feedback from stars 
and UV background, especially the negative feedback ef- 
fects on the SFH and the cosmic reionization history. The 
radiative feedback from PopIII stars includes the ISM pho- 
toevaporation and IGM reionization. The outline of this pa- 
per is as follows. In section 2, we describe the evolution of 
dark matter halos in the ACDM model. Within this frame- 
work, the SFR in a halo can be expressed as an analytic 
formula with some free parameters. Mechanical feedback 
from stars is studied in section 3. In section 4, we intro- 
duce an analytic model to deal with the radiative feedback 
from PopIII. Finally, our discussion and conclusions are pre- 
sented in the last section. Throughout this paper, we adopt 
the cosmological parameters consistent with the 5 years 
WMAP data:»,^/i^ = 0.1358, »a = 0-7 26, f^b = 0.0456, 
h = 0.7, as = 0.812 (|Komatsu et al. 2001 ). 



2. Cold dark matter model and star formation in 
galaxies 

2.1. Redshift evolution of dark matter halos 

In the hierarchical clustering scenario of CDM halos, a 
simple and successful model for the formation and distri- 
bution of spheroidal or ellipsoidal collapsing objects was 
presented in last century based on the theory of Gaussian 
random fields about cosmological density perturbation 
(jPress fc Schechter 1971 iBardeen et al. 19861 ). In this sce- 
nario, gravity governs almost the whole process. The dark 
matter halos increase their own mass throu gh accretin g 
matter and merging with each other lfCacev fc Cole 1993f) . 
Based on this model, Sasaki (1994) proposed an analytic 
formula to describe the formation and evolution of the 
CDM halos. In this formalism, the number density of col- 
lapsed objects with mass in the range (M, M + dM), which 
are formed at the redshift interval (zc, Zc + dzc) and su rvive 
till redshift z, is (jSasaki 1994l : IChiu fc Ostriker 2000[ ). 



N{M, z, Zc)dMdze = aNMizc) 



D(z,)a{M) J D(z) 

dZr 



i/(Zc)(l + Ze) 



dM. (1) 



Here the overdot represents the time derivative, and 
Nm{zc) dM is the number of collapsed objects per unit 
comoving volume within a mass r ange (M, M 4- dM) at 
redshift Zc (jPress fc Schechter 1974 '). In Eq. (P), (5c is a con- 
stant, usually taken to be 1.686 in a matter-dominated flat 
universe (r^m = !)• Thi s value is quite i nsensitive to the cos- 
mological parameters ()Eke et al. 19961 ) . where H{z) is the 



Hubble parameter, -D(z) the growth factor for linear per- 
turbations, and a{M) the rms mass fluctuation on a mass 
scale M . In addition, N{M, z, Zc) in Eq. ([T]) represents the 
formation rate of halos weighted by their survival prob- 
ability. For more accuracy, here we use Sheth & Tormen 

(1999) 's modified formula for the expression of NM{zc)dM , 
which fits the numerical simulations better than the origi- 
nal one, especially at high redshift. Numerical simulations 
also indicate a = 0.707 (see Sheth, Mo & Tormen 2001 for 
more details). 

2.2. Star-formation rate in galaxies and the cosmic SFR 
density 

Stars or even galaxies are confirmed to form through 
coolin g and condensation of baryons within the DM 
halos ( White fc Rees 1978| ). The formation and evolution 
of galaxies and the associated star formation histories 
have been studied extensively via both numerical simu- 
la tions and semi-ana lytic models (Ccn fc Ostriker 1992; 
Chiu fc Ostriker 20001: [Choud hurv fc Srianand 2002 ; 

Springel fc Hernguist 20"ol iNagamine et al. 20061 ). 

Following Gen fc Ostriker (1992), Chiu fc Ostriker 

(2000) and Choudhury fc Srianand (2002), we assum43 
that the SFR in a halo with mass M at z that has collapsed 
at an earlier redshift Zc, is given by 



MsF(M,z,Ze) = /.M^ X iifl^^filexp 



t{z) - t{z,) 



(2) 



Here, 



/*, t{z), tdyn are the fraction of total baryonic 
mass in a halo that will be converted to stars, the 
age of the Universe at redshift z, and the dynamical 
time scale, respectively. In our model, we take fl^^ = 
0.07 (|Choudhurv fc Ferrara 2006') for PopIII stars and = 
0.3 for PopII/I stars. The dynamical time scale ^dyn is given 
by (Chiu fc Ostriker 2000), 



tdyn{z) = 

and here, 

Pvh-{z) = 

Ac(z) = 
d{z) = 

Pc{z) = 



37r 



32Gpvir(-z) 



(3) 



/\c{z)Pc{z), 

187r2 + 82d(z) - 39d2(z), 
n„^{l + zf 



Vlm{l + zf + VlK 

3H^{z) 
8ttG ■ 



1, 



The duration of star formation activity in a halo depends 
on the value of k. Note that k ^ corresponds to the 
star formation occurring in a single burst. By using a 3D 
cosmological code that includes most of the needed physics 
to simulate the formation and evolution of the first galax- 
ies, Ricotti et al. (2002a;2002b;2008) show that the first 
luminous objects ("small-halo objects") are characterized 
by "bursting" star formati on. In this paper, we thus set 
Km = 0.01 for PopIII S F(Samui et al. 20 071) and ku = 
1 for PopII/IdCen fc Ostrike r 1992: Chiu fc Ostriker 20001: 
iChoudhurv fc Srianand 200^) . 



In fact, this is not the most fundamental assumption(see 
Appendix A for details). 
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Consequently, the cosmic SFR per unit comoving vol- 
ume (SFR density, hereafter SFRD) at a redshift z can be 
written as 



Table 1. Model parameters 



Psf(^) 



N{M', z, zc) X Msf(M', z, Ze)dM', 



Ml, 



(4) 

where the lower mass cutoff Afiow(^) at a given epoch de- 
pends on the cooling efficiency of the gas and the different 
feedback processes. Following Barkana & Loeb (2005), the 
lower mass cutoff can be expressed as 



Miow(z) = 9.2x lO^Mg 



16.5km s 



1 



10 



(5) 



where the circular velocity Vc = (2fcBTvir//^"ip)^/^, the 
mean molecular weight /i = 0.61, and Tvir is the virial tem- 
perature of halos. Halos with T^{^ > W^K are able to cool 
via atomic transitions. It is usually considered as standard 
in most semi-analytic models that the minimum mass of 
star- forming halos is Miawiz = 6) ~ 10® M©. When the H2 
content of the gas is considerable, molecular line cooling can 
make gas condense within the small halos (Tyir ~ lO'^K) 
and eventually form stars (iTegmark et al. 19971 ; 

Haiman . Abel fc Rees 20001: iBarkana k Loeb 20051: 

Schneid er" 3t al. 2006t ISamui et al. 20071 ). In the ACDM 
model, the number of halos with masses 10^ — 10* Af© 
(Tvir < lO^K) that are expected to possess PopIII stars at 
high-redshift decr eases after z ~ 10 due to the DM halos 
mergmg process ()Mo fc White 20021) . so these small-halo 
objects can dominate t he galaxy mass fun ction until at 
least redshift z ~ 10 (jRicotti et al. 2002bD . Finally, we 
consider the atomic cooling model (Tvir lO'^K) for 
PopII/I and molecular cooling model (Tvir ~ lO'^K) for 
PopIII. It is notable that the Miow mentioned above 
corresponds to the neutral regions only. 

On the other hand, the Lyman-Werner (912- 
1108 A) backgroimd can also inh ib it the star forma- 
tion ([H aiman, Roc s &: Loeb 19971 : IRicotti et al. 2002al: 
lYoshida et al. 200311 . However, it is interesting that the 
positive feedback of H2 re-formation (e.g. in front of 
HII regions and inside relic HII regions) can counter- 
balance the negat i ve feedback from LW background 
(IRicotti et al. 20011: IRicotti et al. 2002a[ ). As a conse- 
quence, the formation of small-mass galaxies is not 
suppressed, so in this paper we do not consider the LW 
background feedback by assuming that positive feedback 
dominates and star formation in small galaxies is not 
suppressed. 



Parameters Referrence Values 



,11 

^SN 

fill 

J CSC 
rll 

J CSC 



0.05 

0.02 

0.8 

0.1 

0.07 

0.3 



9 - 10 (IRicotti et al . 2002b': 'Sal vaterra fc Ferrara 200l 
IMatsumoto et al. 2005; .Choudhurv fc Ferrara 20061 ). 

Based on the above discussions, we assume a PopIII 
mass fraction, T'ni(z), which means the mass fraction of 
objects forming from gas with Z < Zcrit at a redshift z, i.e . 
the sites of PopIII star formation (jScannapieco et al. 2003h . 
Here T'ni(z) is merely a function of redshift z. The 
transition from PopIII to PopII/I should not happen 
suddenly at a certain redshift, because chemical feed- 
back is a local process: with regions close to star for- 
mation sites rapidly becoming metal-polluted and over- 
shooting Zcnt, and others remaining essentially metal- 
free. PopIII and PopII star forma tion modes could have 
been coeval ( ICiardi fc Ferrara 2005^ ■ so instead of assigning 
T'ni(z) a Heaviside function 0(z — Ztran), here we introduce 
a function: Fiu{f3, z) = {z / Zt^e^nY + [z / ZtrwY], in which 
/3 is a free parameter describing the transition speed. The 
/3 00 limit reproduces the sharp transition case. In this 
work, we take the transition redshift to be ztran — 14, and 
the free parameter (3 — 12, which ensures the mass fraction 
of PopIII is lower than 5 percent at Zond- Now the SFRD 
can be rewritten as 



00 00 



Piii(^) = y dzc j AM'N{M',z,Zc) 

^ Mow 

xMsF(Af',z,Zc) X Tiii(zc) 
for PopIII, and 



(6) 



00 00 



/jii(z)= j dZe j dM'N{M',Z,Zc) 

xMsf(M', z, Zc) X (1 - Tiii(zc)) 
for PopII/I. 



(7) 



2.3. PopIII mass fraction 

With the evolution of the Universe, the character- 
istics of star formation changes from a metal-free, 
massive-star-dominated (PopIII) mode to a metal-enriched, 
normal-star-dominated (PopII/I) mode. Numerical sim- 
ulations indicate that this transition occurs when 
the metal has bee n enriched to a critical value , 
Zcrit ^ IQ-^-^^-'^Zr:, dSromm. Yoshida fc Hernquist 20031: 



ISmith fc Sigurdsson 20071 ). If most of the first gener 
ation stars die as PISNe, the volume-averaged mean 
metallicity would reac h Z ~ 10~'^Z(7) at a red- 
shift of - 12 - 16 (|Yoshida et al. 20041 ). It is sug- 
gested that the PopIII stars are terminated at Zcnd ~ 



3. Mechanical feedback processes 

3.1. Popll/I SNe feedback 

Mechanical feedback is associated with mechanical 
energy injection from SNe explosions and galactic 
winds. Most of works concentrate on the effects of 
the first-generation SNe explosions at very high red- 
shift rather than the win ds from metal-free stars 
(iMori. Ferrara fc Madau 20021 : 



Yoshida et al. 2003 



'Salvaterra et al. 2003; iKitayama fc Yoshida 2005 

Grcif et al, 2007). A consequence of SNe explosions is 
to expel the gas out of the host halo partially (blowout) 
or totally (blowaway) and reduce or empty the reservoir 
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for subsequent star formation. Some numerical simulations 
show that one PopIII PISN = 200Mq (£;sn ~ lO^^gj-gs) 
can deplete its host halo easily(discussed in the next sub- 
section). But there are at least two reasons for an SN 
explosion of PopII/I not being capable of blowing the 
most gas out of its host halo: (i) its explosion energy is 
much lower than that of a PopIII star, and the disperse 
distribution of the lifetime of PopII/I stars cannot ensure 
enough number of the SNe per unit time to deplete the 
halo; (a) a halo with virial temperature Tvir > 10*K in 
the PopII/I epoch(z < 14) having mass of Afvir — lO^Af0, 
corresponding to a binding energy — 10^^ ~ lO^^ergs, 
is able to prevent an SN explosion of PopII/I (gener- 
ally i?sN < lO^^ergs) from expelling most gas out of 
the halo, because for halos with mass Afvir > IO^A/q 
the star formation will not be quenched even if iJgN 
exceeds th e binding energy of halos by 2 orders of 
magnitude ([Kitavama k Yoshida 20051) . Furthermore, the 
enriched metallicity can enhance the gas cooling, which 
can also help to hold the escaping gas back. 

However, the huge energy generated by a (PopII/I) SN 
explosion can, at least, partially heat the cold gas even in 
a high-metallicity environment. This normal SNe feedback 
has been extensively studied both in terms of SNe explo- 
sions and galactic outflows. Due to the feedback of super- 
novae explosions, the gas will be rem oved from the cold 
phase at the rate (jGranato et al. 200^ : 



Afif = -AfsF(Af,z,Zc)e'sN 



^ 500km s 



(7^ 

-IN 2 



(8) 



where tj^^^ ~ 6 x 10 '^Mq ^ is the number of SNe expected 
per solar mass of formed stars, El^^ is the kinetic energy 
of the ejecta from each PopII/I supernova (6 x 10^° erg; 
e.g. Nadyozhin 2003), and eg^^ the fraction of this energy 
that is used to heat the cold gas. Here t]^j^ is evaluated by 
"Starburst99" and adopting a minimum progenitor mass 
of 8A^0 and the Salpeter IMF. Some analyses show that 
above 90 % of the SN kinetic energy may be l ost by radiative 
cooling (jThornton et al. 19981 : [Heckman efal. 20001 . Here 
we set e 



II 

SN 



0.02. As the line-of-sight ve locity dispersio n 
a, we adopt the relationship a ~ 0.65T4 (|Ferrarese 20021 ). 
Finally, we can obtain the SNe feedback on PopII/I SFRD: 



•FB 



(z) = / dz, 



Ml, 



1 - Fm{z,) 



dM' N{M', z, Zc) 



X Mip{M',z,z,). 



(9) 



where Afgp(Af, z, Zc) is the regulated SFR by the feed- 
back from SNe explosions in PopII /I dominated halos (with 
Tvir > lO^fsT, see Appendix A for detail). 

3.2. Mechanical feedback from first-generation stars 

The ultima te fate of a metal- free star depends criti cally 
on its mass (|Heger fc Wooslev 20021 : iHeger et al. 20031 ): (1) 
8Mq < Af* < 25Af0 (these stars explode as core-collapse 
SNe and leave neutron stars behind), (2) 25Af0 < Af* < 
4OAf0 (these explode as faint Type II SNe and leave black 



holes behind), (3) AOMq < A'h < UOMq (these do not 
explode as SNe and directly collapse into black hole^ , (4) 
14OAf0 < A/, < 26OAf0 (these explode as PISNe, causing 
complete disruption), (5) 26OAf0 < Af, (these collapse, in 
the absence of rotation, directly into black holes), so not 
all the PopIII stars can die as an SN explosion. We employ 
a slightly top-heavy IMF (Larson, 1998) for PopIII stars 
(5OM0 < A'U < 500Mq): 



dN 



(dlogAf*; 



cx (l-f Af^/Afc)- 



-1.35 



(10) 



where Afc — 100 A/0 is the characteristic stellar mass of 
PopIII. As a result, the PISN is almost the only type of su- 
pernova explosion for PopIII. Bromm, Yoshida & Hernquist 
(2003) show that, for a halo of mass M ~ 10^ Mq at z 20, 
a PISN of mass Af, = 25OAf0(£'sN ~ lO^^ergs) can disrupt 
the halo completely. A similar result has been obtained 
by Grcif et al. (2007), who find that a PISN with mass 
Af, — 2OOAf0(i?gN ~ lO^^ergs) can deplete the whole host 
halo. 

The photoevaporation effe ct might be particularly im- 
portant for PopIII objects (jCiardi fc Ferrara 20051 ). In 
small-halo objects(Tvir < W^K), photoevaporation alone 
from OB stars can produce strong galactic winds that ex- 
pel most of the gas from galaxies. Galactic winds pro- 
duced by an SN explo sion may be important after about 
10 Mvr(Ric otti et al. 20081. Because the internal photoe- 
vaporation in a PopIII object is able to deplete most of the 
gas before the first SNe explode, in fact, the first SNe explo- 
sions have no chance to exert negative feedback on the star 
formation in these halos. Finally we only consider the SNe 
feedback in halos of Tvir > W^K in this paper. For those big 
halos(Tvir > 10"^ K), the high coohng effect and deeply grav- 
itational well can confine the gas photoevaporation. And 
even the most powerful PISN ( E'sN — lO^^erqs) cannot sig - 
nificantly blow the gas awav l|Kitavama fc Yoshida 2005f ). 
In reference to the mechanism of the SNe feedback of 
PopII/I, we use the same formula but with different pa- 
rameters to describe the PopIII-SNe feedback in these halos 
{M > IO^Mq), 



^iii piii 

TiVSN _ ^ aV (A f -y -y will %N^SN 
^^^iil — o^^^SFiAi, Z, Zcl ■ 



4.5e™AfsF(Af,z,Zc) 



^ 500km s 



-IN 2 



(11) 



Here 



1.5xlO~^Ar0^ (evaluated with 
the PopIII IMF and mass range). For lack of information 
about the strengt h factor of the SNe feedback of Pop III, 
we set e^'i^ = 0.05 (|Granato et al. 200l:lLai3i et al. 2006f ) as 
our reference value. One can obtain the PopIII SFRD with 
SNe feedback: 

Pmiz) ^ J dzcFniizc) J N{M',z,Zc) 

xMl^^{M',z,Zc)dM' (12) 

where A^gp(Af, z, Zc) is the regulated SFR by the feedback 
from SNe explosions in PopIII halos of Tvir > 10"'-fi'(see 



^ But some of them experience a pulsating instability and 
eject their outer envelope, again leaving black holes behind. 
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Fig. 1. Evolution of SFRD with the redshift z. Left panel: the crossed data points are taken from Hopkins & Beacom 
(2006). The squared ones come from Reddy et al. (2008). Three circled points are from the UDF data (Bouwens et al. 
2005). The dashed line represents the PopII/I SFRD without SNe feedback, while the dot-dashed one with the SNe 
feedback. Right panel: the original SFRD of PopIII (Tvir > lO'^K, the dot-dashed hnc) is suppressed down about one 
order of magnitude by the SNe feedback (the dotted line). 



Appendix A for detail). Of course, it is definite that the 
number density of halos with Tvir > W^K is tiny in PopIII 
epoch according to the hierarchical clustering scenario. 

Figure [T] shows the SNe feedback in PopIII (Tvir > 
lO^K) and PopII/I halos. In the left panel, we scale down 
all the observation data by a factor 1.9 to make it con- 
sistent with the IMF used in our model (Mf '"^O.SM©, 
Mf =^'==100Mo for PopII/I stars and Salpeter type IMF). In 
the right panel, the original SFRD of PopIII(Tvir > lO^K) 
is cut down about one order of magnitude by the SNe feed- 
back. 

Figure [Umanifests the relative feedback strength of SNe 
feedback[pfjf(z)/piii(2;) or pjy^{z) / pii/i{z)] reduces as the 
redshift decreases, because massive halos appear abundant 
only at low redshift according to the "bottom-up" hierarchi- 
cal structure model (see Mo & White 2002 for details) and 
massive halos whose binding energy is much higher than 
before can reduce the feedback effect from SNe explosions. 
On the other hand, massive stars disappear as the metal- 
licity of gas is enriched. Small star formation with a long 
lifetime and small £'sn is enhanced in the metal-enriched 
cloud. 

4. Radiative feedback and cosmic reionization 

The ionizing radiation produced by massive stars or quasars 
can have local effects or long-range effects, either affecting 
the formation and evolution of nearby objects or joining 
the radiation produced by other galaxies to form a back- 
ground. Mini-QSOs that are usually thought to be powered 
by accretion onto the first black holes may provide an X- 
ray background when the universe has a certain number of 



mini-QSOs in a unit comoving volume. This background 
would ionize the HI to iJ+ and by direct photoion- 
ization or coUisional excitation from a high-energy pho- 
toelectron produced by photoionization. If this event oc- 
curs at the core of a protogala xy, the fraction of molecular 
hydrogen would be p romoted (|Haiman. Abel fc Rees 20001 : 
iGlover fc Brand 20031) via 

H + c- ^ H- -^7, (13) 
H- + H ^ H2 +e-. (14) 

Dijkstra et al.(2004) and Salvaterra et al.(2005) show that 
the hard X-ray from the same sources will produce a 
present-day soft X-ray background. They find that the 
models with accreting black holes will overproduce the ob- 
served X-ray background by a large factor. A population 
dominated by mini-QSOs could still partially ionize the 
IGM at z > 6, but its contribution can be severely con- 
strained if the X-ray background is resolved further into 
discrete sources. By considering the constraint from the 
soft X-ray background, Ricotti & Ostriker (2004b) show 
that mini-quasars at high-z do not overproduce the X-ray 
background and can still produce a significant contribution 
to reionization and the optical depth of electrons t^. 

Based on the investigation of the environment and 
reionization process around the highest redshift QSOs hav- 
ing Gunn-Peterson troughs {z > 6.1), Yu & Lu(2005) ar- 
gue that a significant fraction of hydrogen in the Stromgren 
sphere around QSOs is ionized by photons from stars and 
that only about several percent to at most 10% — 20% of 
the total hydrogen is left (e.g., in minihalos, halos, or high- 
density subregions) to be ionized by QSO photons. Willott 
et al.(2005) show that the current constraints on the quasar 
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population give an ionizing photon density <C 30% that of 
the star-forming galaxy population by analyzing the obser- 
vational data from the Canada-France-Hawaii Telescope. 
They conclude that active galactic nuclei make a negligible 
contribution to the reionization of hydrogen at z ~ 6. But 
this argument only applies to quasars at z < 6, one can 
still have a high— z population of mini-quasars that par- 
tially ionize the I GM at z > 10 witho ut overproducing the 
X-ray backgroun d (|Ricotti et al. 20051 ). 

Pre-ionization by X-rays can increase the IGM opti- 
cal depth from Te — 0.06 given by stellar sources only 
to 0.1 < Tc < 0.2 (|Ricotti fc Ostriker 2 004b). From 4-5 
years ago, it was necessary to fit an observational value of 
To ~ 0.17 measured by WMAP satelhte. Recent WMAP 
data indicating 0.068 < Te < 0.1 maybe imply that pre- 
ionization by X-rays is not as important as before. In this 
paper, we focus on the ionization by stellar sources and ig- 
nore the effects from QSOs or mini-QSOs. We will consider 
the AGN feedback carefully in the work in preparation. 

4.1. Photoevaporation in small-halo objects (T^n- < IQ'^K) 

PopIII stars quite likely reside in small DM halos and a 
metal- free environment. It benefits the photons escaping 
and the radiative cooling. As mentioned in the last sec- 
tion, photoevaporation is important in small-halo objects 
(Tvir < 10^ K, e.g. most of the PopIII halos). It can prevent 
the size of HII regions from exceeding R^if — 5/i~^kpc, 
the mean free path of ionizing photons, about the size 
of the dense filamen ts and the virial radii of the halos 
(jRicotti et al. 2002bl ). When the HII regions become bigger 
than the filaments, molecular hydrogen is destroyed and the 
star formation is suppressed. An analytic model has been 
developed t o describe the propagat ion of ionization fronts 
in the IGM (|Barkana fc Loeb 200lh : 



< n 



ISM 



> 



dV„ dN^ 



dt 



dt 



Q!B 



(15) 



where is the number of ionizing photons produced by 
the source. Here Vp is the ionized proper volume, ae — 
2.6x 10~^'^c m'^ s~^ is the case B recombination coeffici ent at 
T ~ lO'' K (jSeager et al. 19991 : iBarkana fc Loeb 200lh . and 
'^ISM ^® proper number density of the hydrogen atoms. 
We assume the ISM in halos is homogeneous, which means 
< («ism)^ >=< (j^ism) so the analytical solution is 



1 



dN. 



7 



X exp < — ttB 



wher^l 



< <SMit') > dt' 

<7i|m(^") >dt"]dt' 



< nfsuit) >- A,(i) X ^ X 4 X (1 + z{t)f, 

0:022) ' 



(16) 



(17) 



n° =1. 



X 10" 



and is the present number density of hydrogen. As men- 
tioned in Eq. ([3]), Ac is the overdensity relative to the crit- 
ical density at the collapse redshift. Following Barkana fc 



^ The mass fraction of baryons in a halo is , so the baryon 
density in halos is at least Ac{t) x times the one in IGM. 
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Fig. 2. The duration time of star formation tfb = t' ~t'{zc) 
versus the masses of halos. 



Loeb (2001), we can evaluate dN^/dt by using 
,111 



dN.y 
"df 



a - 1 
a t.. 



(18) 



where a is the index in the mass-luminosity relation (for 
OB stars a = 4.5), and n^^ the number of ionizing pho- 
tons released per bary on of stars formed ijSchaerer 20031 : 
iHaiman fc Brvan 2006f ). Fang fc Cen (2004) present a rela- 
tionship between n^^ and the stellar mass, M^,. Because 
is not sensitive to M* in the mass rang e of concern, we se t 
as a fixed value 80000 in this paper ljFang fc Cen 2004[ ). 
Most massive stars fade away with the characteristic tim e 
scale ts = 3x 10^vr (|Bond et al. 1984 lAlvarez et al. 20061 ). 
For a halo of mass M, star formation initiated at t{zc), the 
proper size(or radii) of the HII region at later time t, is 



i?(M,t(z.),i) = ^ A/m^miMi 

47r S2n,771r, 



X Vp{t,t{z,)) 



1/3 



where flU is the escape fraction of ionizing photons from 
the sources(or resolution elements in 3D numerical simula- 
tions, see Ricotti et al. 2002a, 2002b, 2008), and /™ ~ 
0.3 ~ 0.9 is suggested by Choudhury fc Ferrara (2006) 
and Mao et al.(2007). Yoshida et al.(2007) provide a time- 
averaged ionizing photon escape fraction as a function of 
stellar mass M*. For the mass range adopted in this work, 
we use a fixed value flU = 0.8. As mentioned above, star 
formation will be suppressed when R{M,t' (zc),t') is equal 
to the size of the filaments. Thus the relation between 
t' ~ t'{zc) and M is given by 



R{M,t'{z,),t') 



i?Hn/(l + ^')- 



(20) 



We illustrate this relationship in Fig. [21 where tfb — t' — 
t'{zc) is the duration time of star formation. We plot five 
cases Zc = 10, 15, 20, 25, 30, which are actually not straight 
lines in the figure. For the halos with the same mass, tfb 
increases as the redshift decreases. The bigger the halo, the 
less time it spends enlarging the HII region to the size of 
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Fig. 3. Evolution of SFRD as a function of the redshift 
z with/ without the photoevaporation feedback effect. The 
sohd hne represents the SFRD without feedback. The dot- 
ted one stands for the case that the SF is suppressed by pho- 
toevaporation. Both of them represent the PopIII SFRD in 
halos with Tvir < W^K. 



the filaments. If tfb < tg, the photoevaporation feedback 
occurs before the first SN exploded. And the inequality 
c X Tfb > -RhhV (1 + the speed of hght) is valid for all 
these cases, which means the speed of the propagation of 
the ionization fronts cannot exceed the speed of the light|f|. 

Now we can calculate the SFRD in PopIII halos with 
Tvir < lO^K by using 



oo Af(zc40 AT) 

pfif(z,< lO^i^) = ydzc^^iii(zc) J N{M',z,Zc) 

z M{zc,10^K) 
xe{t{z,) + TFB - tz) X Msf{M', Z, Ze)dM'(21) 

where 6* is a Heaviside function. Figure [3] shows the ef- 
fects of this kind of feedback: (i) they suppress the PopIII 
SFRD down about one order of magnitude at the peak 
(z ~ 15);(m) photoevaporation gets stronger when redshift 
decreases. And we also find that the feedback will become 
weak when fill is tuned down. Moreover, the SFRD is not 

TTT 

sensitive to the /^^^. 

4.2. Radiative feedback in large halos (T^i,. > W^K) 

For halos with Tvir higher than the cooling temperature 
of the hydrogen atoms(10'*iir), gas cools via HI emission 
lines, and the photoevaporation by internal sources cannot 



However, it seems that the offcenter sources may make a 
superluminal spread when their HII regions get an overlap. But 
it is not real superluminal behavior. 



expel gas significantly. But the ionizing background can 
still suppress some small halos in which the star forma- 
tion is about to initiate. Due to the IGM reionized by UV 
photon escape from star-forming galaxies, the temperature 
of the gas in halos is enhanced. This can dramatically in- 
crease the Jeans mass. Furthermore, numerical simulations 
indicate that the photoionizing background can completely 
suppress galaxy formation in halos with circular velocity 
Vc < 35 km s~^, while the mass of cooled baryons is re- 
duced by 50% for halos wi th circular velocities ^ 50 km s~^ 
(jThoul &: Weinberg 19961 ). In the ionized fraction of the 
universe, we assume complete suppression of star forma- 
tion in halos below circular velocity = 35 km s~^ and 
no suppression above circular velocity of 95 km s^^. For 
intermediate masses, we adopt a linear fit from 1 to for 
the suppression factor (as in Bromm & Loeb, 2002). Thus 
in the HII region the SFRD can be expressed by 



dzc J dM' N{M',z,Zc) 



1 - Fulizc) 



X Mip{M',z,z,) X WiVo). 
for PopII/I and 

oo oo 

> lO^X) = y" dzcFiii(zc) J N{M',z,z,) 



(22) 



xM^P(M',0,Zc) X M^(Vc)dM(23) 



for PopIII, where VF(14) is 
'0, 

-35)/(95 - 35), 




Vc < 35km s^i 
35 < K < 95km s^i (24) 
Vc > 95km s^i 

We plot two extreme cases in Fig. |4| the stars form in 
the completely reionized universe in HII region and in the 
completely neutral universe in HI region. This radiative 
feedback from Popll/l{left panel) and PopIII has the same 
behavior: feedback gets stronger when redshift increases. 
Similar to the SNe mechanical feedback mentioned above, 
the more massive the halo, the less the SF can be affected. 

For the ionized fraction of the universe at a given red- 
shift z, it can be described by the following equation. For 
simplicity, we assume that all the Lymann continuum pho- 
tons escape from a star-forming galaxy are involved in 
reionizing the IGM, and the photons emitted from the 
sources will immediately join the action with the atom, re- 
gardless of the photon propagation between the source and 
the atorr[f|. Then the fraction of io nized hydrogen Qmi, th e 
so-called filling factor, evolves as (jBarkana &: Loeb 200l[ ). 



dg 



HII 



dz 



{l + zfnl ^ d-z - + "H^™^ dl' 

(25) 

where the volume-averaged clumping factor of the IGM, C, 
is defined as C = (n^) / ((1 -I- z)'^ny)^. The first term on the 
right is the rate of ionization and the second term is the 



Recent work shows that the finite speed of light will lead to 
a substantial change of the growth rate in the ionized volume 
(see Qiu et al. 2008 for more details). 
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Fig. 4. Evolution of SFRD in two extreme cases: the stars form in the completely reionized universe (HII region, the 
dashed line) and in the completely neutral universe (HI region, the dot-dashed line in the lefthand panel and the dotted 
line in the right one). Observational data are the same as Fig.[TJ 



rate of recombination, weighted by the Qhii, as recombina- 
tions take place only in the ionized region. Then dN-y/dt is 
obtained from the SFRD calculation as, 

dt f*mp ^' 

where psF{z,T^ir > lO^K) is the SFRD in halos of T^i, > 
1Q*K (including all PopII/I halos and a few PopIII ha- 
los) at z. Because the strong photoevaporation feedback in 
small-halo objects can quench the SF, which stops the ion- 
izing photons to escape to the IGM. This part of SFRD 
contribute little to the IGM reionization. The value of rij 
depends on the IMF of the forming stars. For a Salpeter 
IMF(with 0.5 Mq < < 100 Mq), nl^ is about 5,400, 
evaluated by "Starburst99" with metallicity Z = 0.008, 
however, for the metal-free stars, the IMF could be biased 
towards very massive stars. As mentioned in the last sec- 
tion, we still used n}^^ ^ 80, 000 as our reference value. For 
the clumping factor of IGM, C, we adopted the simple form 
given by Haiman & Bryan (2006). 

By solving Eq. (P5|) (see Appendix B for detail), one 
can find that SF and reionization affect each other. When 
star formation gets high, HI ionization becomes strong 
and star formation is suppressed. Hopefully, JWST can 
help us to understand their relationship (SF and reioniza- 
tion) exactly. The result is plotted in Fig.EKZe/i panel). 
One may find that the solid line separates from the 
SFRD(PopII/I-|-PopIII) in HI region and eventually merges 
into the SFRD line in HII region as the filling factor Qhii 
(see Fig. ^ increases from zero to unit. The PopIII SFRD 

^ http:/ /www. stsci.edu/science/starburst99/ 



in halos of Tvir < W^K is taken into account in the 
righthand panel in Fig. O By comparing these two panel, 
we can find that the main part of PopIII SFRD is still 
contributed by those stars in small-halo objects. Since the 
PopIII PISNe contaminate the IGM with heavy elements 
at very high efficiency (recall the ultimate fate of metal- 
free stars), PopIII ends its era by itself quickly (survival 
time ~ 2 X lO^yr.). The ignition of PopII stars makes the 
reionization continue and eventually completes it at z ^ 6. 

The reionization history and the optical depth of elec- 
trons are plotted in Fig. [S] Here we take the singly ionized 
Hell into account in righthand panel in Fig. [6l The equa- 
tion for calculating the optical depth of electrons is: 

f ds 

Tc{z) = / nc{z)ac—dz (27) 
J dz 

where n^ (z) is the mean electron number density in the uni- 
verse, fjc stands for the Thompson cross section for electron 
scattering, and ds — ~cdz/H[z){l + z) is the proper line el- 
ement. Figure iniillustrates that the model output can fit the 
observation from the highest QSOs (zre ~ 6), and WMAP 
5 years resuh (to = 0.0841°;°}^ Komatsu et al. 2008) at the 
same time. Most of the UV photons emitted from small ha- 
los that have the overwhelming majority in PopIII objects 
are confined in As shown in the lefthand panel in 

Fig. [6l PopIII stars are not able to reionize the universe sig- 
nificantly (about 20% at z ~ 14). In the righthand panel, 
the optical depth to Thomson scattering provided by PopIII 
stars is less than 22%. 

5. Discussion and conclusions 

In this paper we develop an analytic model of the mechani- 
cal and radiative feedback of star formation based on some 




Fig. 5. Evolution of SFRD with feedback. Right panel: the solid line stands for the SFRD in halos of Tyir > W^K. Other 
lines are the same ones as in Fig. 2] Left panel: The PopIII SFRD in halos of Tvir < lO'^K. 




Fig. 6. Evolution of the filling factor Qhii {left panel) and the optical depth to Thomson electron scattering Te {right 
panel) with the redshift. 



results of numerical simulations. Basically, our results can 
fit the observational data. However, this model also de- 
pends on some initial parameters (most of them are listed 
in Table 1) and working assumptions. 

(i) It is possible to fit the reionization result from the 
highest QSOs observation by providing another pair of fl^ 
and fll^. This is because the IGM reionization mainly de- 
pends on nl^ and fl^ x fH^. The value of n}^ is nearly fixed 
under certain conditions (e.g. metallicity, IMF, and stellar 
mass range), but lacking the SFRD observational data con- 



straint on fl \ one can find many pairs of fl and fH^ by 
keeping fl^ x fH^ invariant to guarantee To = 0.084to;g}g. 

(ii) As for the assumption on Fiii{z), there are a lot 
of similar functions. This is an uncertainty factor in this 
model. Hopefully, it may be determined by the numerical 
simulations and even the future JWST observation. 

{Hi) We only consider the negative feedback in this 
work, the neglect of some positive feedback, such as the X- 
ray background or dense shell (from SNe explosions) frag- 
mentation, results in the underestimation of the SFRD, 
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Qhii7 ''e, and so on. This is what we wiU investigate in 
the next work. 

(iv) Although we can fit the optical depth to Thomson 
scattering in this work, it does not mean the X-ray pre- 
ionization is not important. Obviously, our results are close 
to the low fringe of the error bar. Perhaps, the left space 
in the error bar is waiting for the contribution from X-ray 
pre-ionization. 

In the following, we summarize the results of this work 
in more detail. 



1. The SNe feedback can 
tion in large halos (T^ir > 



suppress the 
lO'^i^). And 



star forma- 
the relative 
or 



feedback strength of SNe feedback [pf^ (z) / piu{z) 
Piyii^) / Pii/ii^)] reduces as the redshift decreases, because 
massive halos appear abundant at low redshift and massive 
halos whose binding energy is much higher than before can 
reduce the feedback effect from SNe(see EqEland Eq fTT|) . 

2. The radiative feedback mechanism can affect the SFR 
in both small (Tvir < W'^K) and large(rvir > lO'^K) PopIII 
halos. In small halos, the shallow gravitational well and 
the poor cooling-efficiency cannot prevent the hot gas from 
escaping due to the strong photoevaporation from the mas- 
sive OB stars. The radiative feedback can suppress the 
PopIII SFR considerably. In large halos, the ionizing back- 
ground enhances the Jeans mass and makes the SFR sup- 
pressed more or less. 

3. The SFRD in smaU PopIII halosjTvir < lO^'X) 



sensitive to 
ously. If /, 



TpB, which depends on fUl and 



III 

CSC 

the 



and 



IS 



sen- 



increases, TpB will decreases and 
so does the SFRD. Furthermore, our results support 
this viewpoin t: early star format ion is likely to be self- 
regulated (iRicotti et al. 2002al : iRicotti et al. 2002b : 



Ciardi fc Ferrara 2005 



Yoshida et al. 2003': 

Kitavama & Yoshid a 20051 : IRicotti et al. 20081) . 

4. Although PopIII stars have high value of /™ and 
n}^^, they are not able to reionize the universe considerably 
(about 20% at 2; ~ 14). And the optical depth to Thomson 
scattering provided by PopIII stars is less than 22%. 

In general, we find the radiative feedback is important 
to the formation of the early generation stars. It suppresses 
the star formation considerably. But the mechanical feed- 
back from the SNe explosions is not able to affect the early 
star formation significantly. The radiative and mechanical 
feedback dominates the star formation rate of the PopII/I 
stars. The feedback on SFRD from first generation stars 
is very strong and should not be neglected. However, their 
effect on the cosmic reionization is not significant, which 
results in a small contribution to the optical depth of elec- 
trons Tp. 
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Appendix A: How to get the SFR formula with 
SNe feedback 

Following Cen & Ostriker (1992), a halo with mass M has 
the initial baryonic gas Mb(0) = |p- x M. Later on, partial 
baryonic gas condenses into cold gas Mcoid- It is reasonable 
to assume that the rates of the mass change are propor- 
tional to Mb(t) at time t: 



dMb(0 
d< 



^Mb(t) and 



det 



+Mb(t). (A.l) 



sookms"^ 



; therefore, 



M.'SMi) - Mcoid(0 - f-' I MsF{t')dt' - I A/if (t')clt' 



= Afcold(t) - 







500 \ ■ 



Vc J 



MsF{t')dt', 



where the circular veloc ity is a functio n with 
two variables M and z ([Barkana &: Loeb 20011 ). Then 

1 /3 

{fFtzJTf^} varies between 0.65 and 1.0 when z is be- 
tween 30 and zero. For simplicity we use its intermediate 
value 0.75 in this work: 



V'c(Af,Zc) = %/0?75 X 23.4 



M 



lO^h-^Mr. 



1/3 



1 + Zc 



1/2 

(A-7) 



Solving Eq. (fO)) again with AfsF(O) 
the SFR with PopII/I SNe feedback: 
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0, one can get 



This is also an implicit assumption in Cen & Ostriker 
(1992). By solving Eq. (jA.l[) with the initial condition, we 
obtain 



Mb{t) = A/b(0) exp(-t) and 
Meoid(0 = / AUt')dt' = Mb(0) [1 - exp(-t)] 



(A.2) 



Cold gas can form stars with efhciency /*. And we as- 
sume that the newly formed star mass per unit time, Msf, 
is proportional to the net mass of cold gas at that time 



dt 



exp 




(A., 



where Su{M,Zc) = 1.82 x 10^ e^^^M'^/^{l + z^)-^ for 
the feedback from PopII/I stars, while Sui{M,Zc) — 
7.58 X 103e|^]^M-2/3(l + z^)'^ for the feedback from 
PopIII stars in halos with Tvir > W^K. After substitut- 
ing [t{z) — t{zc)] I [Ktdyn(zc)] for wc havc the SFR with 
SNe feedback: 



MsF(i) = /* X ACi^i), 
here 



Ml 



S*d(t) = Meoid(<) - 



MsF(i')di'. 



(A.3) 



(A.4) 



Solving Eq. jXll with Msf(O) = (because M"oid(0) = 0)^ 
we get 



dMsF(t) 
dt 



/,A/b(0) X texp(-t). 



(A.5) 



If we replace t with \t(z) — t{z^\ j \Kt^y^-J^z^\^ one can 
easily obtain Eq. ([2|). Since the feedback of supernovae 
explosions removes the cold gas at the rat43 Affj^ = 



In halos of mass M > 10 or virial temperature of Tvir > 
10" is: after z = 30, even a PIS N with Esn = lO'^^ergs ca nnot 
make a substantial gas outflowl jKitavama fc Yoshida 20051 ). So 
we assume the baryons loss in a halo should not be taken into 
account. 



M^i^{M,z,z^) = 



llbM 



exp 



- exp 



17m5'll(M, Zc)Ktdyn 

A"5„(Af, Ze) + 1 



i(z)-i(Zc) 



'*^dyn(Zc) 

i{z)~t{z,) \\ 

Ktdyn{Zc) J I 



(A.9) 



and 



Ml]^{M, Z, Ze 



■ exp 



^^mS'lIl(Af, Zc)K.tdyn 

fl''Sm{M,z, 



exp 



t(z) - t(Ze) 



'^idyn(Zc) 

t{z)-t{zc) y 

'^^dyn(Zc) J 



(A.IO) 



Obviously, Eq. (|A.9|) and Eq. (jA.lOp reduces to Eq. ^ 
when S'iii(Af, Zc) approaches to zero. 
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Appendix B: Radiative feedback from massive stars 



We focus on the radiative feedback from the ionizing back- 
ground. Similar to Samui et al. (2007), we describe the 
strength of radiative feedback via the circular velocity 
of halos Vc and the fraction of ionized hydrogen Qhii- 
Considering the IGM in a unit volume after the ignition of 
the first stars, partial volume was reionized by the UV pho- 
tons from massive stars in the massive halos(Tvir > 10^ K). 
For the ionized part of the halo, their SFR will be sup- 
pressed by this ionizing background, which also depends 
on their mass. On the other hand, the SF will continue in 
the rest of the volume, where the IGM is neutral, so the 
total SFRD is contributed by these two parts at the same 
time: 



Psf(^, > lO^K) = (z) (1 - Qhii(^)) + pfp (^)Qhii(^) 
-pWiiz,> 10^ K)\il~Qmi{z)) 



Pn"W+Pm(^,>10"if) 



•HII 



>HIl(^)- 

(B.l) 



where pf/(z) stands for pfi^(z), pg{(z) for pf^ j (z, > W^K). 

By combining Eq. Eq. and Eq. (|b.1|) . one can 
get, 



dQmijz) 
dz 

here, 



(B.2) 



- ± J „ni .III ^dp™(z, > lO^i^) dp™(z, > lO^K) 

^ ' ~ 7^0 ] 'S •'esc 



dz 



dz 



,11.11 fdpW'iz) dpf/(z) 



7 J CSC 



dz 



dz 



aB{l + zfn'HCiz)^, 
dz 



„i,„idpHl(z,>104i^) , dp™(z) 



J ( 



■'7 J CSC 



dz 



7 •/ CSC 



dz 



